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Abstract 



By combining (i) the economic theory of rational expectation bubbles, (ii) behavioral finance on imitation and herding 
of investors and traders and (iii) the mathematical and statistical physics of bifurcations and phase transitions, the 
log-periodic power law (LPPL) model has been developed as a flexible tool to detect bubbles. The LPPL model 
considers the faster-than-exponential (power law with finite-time singularity) increase in asset prices decorated by 
accelerating oscillations as the main diagnostic of bubbles. It embodies a positive feedback loop of higher return 
anticipations competing with negative feedback spirals of crash expectations. We use the LPPL model in one of its 
incarnations to analyze two bubbles and subsequent market crashes in two important indexes in the Chinese stock 
markets between May 2005 and July 2009. Both the Shanghai Stock Exchange Composite index (US ticker symbol 
SSEC) and Shenzhen Stock Exchange Component index (SZSC) exhibited such behavior in two distinct time periods: 
1) from mid-2005, bursting in October 2007 and 2) from November 2 008, bursting in the beginning of August 2009. 



We successfully predicted time windows for both crashes in advance (iSornette , 



2007 



Bastiaensen et al 



the sam e methods used to successfully predict the peak in mid-2006 of the U S housing bubble (IZhou and Sornette . 



)9) 



2009) with 



2006bl) and the peak in July 2008 of the global oil bubble dSornette et al 



2009). The more recent bubble in the Chinese 



indexes was detected and its end or change of regime was predicted independently by two groups with similar results, 
showing that the model has been well-documented and can be replicated by industrial practitioners. Here we present 
more detailed analysis of the individual Chinese index predictions and of the methods used to make and test them. 
We complement the detection of log-periodic behavior with Lomb spectral analysis of detrended residuals and (H, q)- 
derivative of logarithmic indexes for both bubbles. We perform unit-root tests on the residuals from the log-periodic 
power law model to confirm the Qrnste in-Uhlenbeck pr operty of bounded residuals, in agreement with the consistent 



model of 'explosive' financial bubbles (ILin et al 



2009) 
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1. Conceptual framework and the two Chinese bubbles of 2005-2007 and 2008-2009 



The present paper contributes to the literature on financial bubbles by presenting two case studies and new empir- 
ical tests, in support of the proposal that (i) the presence of a bubble can be diagnosed quantitatively before its demise 
and (ii) the end of the bubble has a degree of predictability. 

The s e two claims are highly contentious and collide against a large consensus both in the academic literature 



dRossei 



bubbles. 



2008 ) and among professionals. For instance, in his recent review of the financial economic literature on 



Gurkavnak 



(120081) reports that "for each paper that finds evidence of bubbles, there is another one that fits 
the data equally well without allowing for a bubble. We are still unable to distinguish bubbles from time-varying or 
regime-switching fundamentals, while many small sample econometrics prob lems of bubble test s remain unresolved." 



Greenspan! (120021) . at a summer confer- 



Similarly, the following statement by former Federal Reserve chairman Alan|j 
ence in August 2002 organized by the Fed to try to understand the cause of the ITC bubble and its subsequent crash 
in 2000 and 2001, summarizes well the state of the art from the point of view of practitioners: "We, at the Federal 
Reserve recognized that, despite our suspicions, it was very difficult to definitively identify a bubble until after the fact, 
that is, when its bursting confirmed its existence. Moreover, it was far from obvious that bubbles, even if identified 
early, could be preempted short of the Central Bank inducing a substantial contraction in economic activity, the very 
outcome we would be seeking to avoid." 

To break this stalemate, one of us (DS) with Anders Johansen from 1995 to 2002, with Wei-Xing Zhou since 2002 
(now Professor at ECUST in Shanghai) and, since 2008, with the FCO group at ETH Zurich (www. er . etttz . ctt/f co/) 
have developed a series of models and techniques at the boundaries between financial economics, behavioral fi- 
nance and statistical physics. Our purpose here is not to summarize the corresponding papers, which explore many 
different options, including rational expectation bubble models with noise traders, agent-based models of herding 
traders with Bayesian updates of their beliefs, m odels with mixtures of nonlinear trend followers and nonlinear 



value investors, and so on (see 



Kaizoii and Sornettd (120091) : 



Sornette (2003b) and references therein until 2002 and the two recent reviews in 



Sornette and Woodardl (|2009) and references therein). In a nutshell, bubbles are identi- 



fied as "super-exponential" price processes, punctuated by bursts of negative feedback spirals of crash expectations. 
These works have been translated into an operational methodology to calibrate price tim e series and diagnose bubbles 



2003bj) and more recently successful 



as they develop. Many cases are reported in Chapter 9 of the book d Sornettd, 
applications have been presented with ex-ante public announcements posted on the scientific international database 
arXiv . org and then published in the referred literature, which include the diagnostic and identification of the peak 
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in mid-2006 (Zhou and Sornett 



Kindlebergei 



10) 



2000) and 



2006b|). and the oil price peak in July 2008 dSornette et al 



time of the b ubble for the UK real-estate bubble in mid-2004 ( IZhou and Sornette , 



2003a), the U.S. real-estate bubble 



2009). 



Sorne tte (2003b|) have identified the following generic scenario developing in five acts, 
which is common to all historical bubbles: displacement, take-off, exuberance, critical stage and crash. For the 
Chinese bubble starting in 2005, the "displacement" and "take-off" can be associated with the split share structure 
reform of listed companies in 2005. Before the reform, only about one third of the shares of any listed company in the 
Chinese stock market were tradable. The other two-third shares were non-tradable (not allowed to be exchanged and to 
circulate between investors), and were owned by the state and by legal entities. The tradable stocks acquired therefore 
a significant liquidity premium, and were valued much higher than their non-tradable siblings, notwithstanding the fact 
that both gave the same privilege to th eir owners in terms of voti ng rights and dividends. Since 2001, the Chinese stock 



market entered an anti -bubble phase dZhou and Sornette , 



2004) with the Shanghai Stock Exchange Composite index 



falling from its then historical high 2245 on 24 June 2001 to the historical low on June 6, 2005. On 29 April 2005, 
the China Securities Regulatory Commission launched the pilot reform of the split share structure. The split share 
structure reform is defined as the process to eliminate the discrepancies in the A-share transfer system via a negotiation 
mechanism to balance the interests of non-tradable shareholders and tradable shareholders. On 4 September 2005, the 
China Securities Regulatory Commission enacted the Administrative Measures on the Split Share Structure Reform 
of Listed Companies^, which took effect immediately. It is widely accepted that the split share structure reform was a 
turning point which triggered and catalyzed the recovery of the Chinese stock market from its previous bearish regime. 
For the Chinese bubble starting in November 2008, the "take-off" can be associated with China's policy reaction on the 
global financial crisis, with a huge RMB 4 trillion stimulus plan and aggressive loan growth by financial institutions. 

Here, we present an ex-post analysis of what we identified earlier in their respective epochs as being two significant 
bubbles developing in the major Chinese stock markets, the first one from 2005 to 2007 and the second one from 2008 
to 2009. The organized stock market in mainland China is composed of two stock exchanges, the Shanghai Stock 
Exchange (SHSE) and the Shenzhen Stock Exchange (SZSE). The most important indices for A-shares in SHSE and 
SZSE are the Shanghai Stock Exchange Composite index (SSEC) and the Shenzhen Stock Exchange Component 
index (SZSC). The SSEC and SZSC indexes have suffered a more than 70% drop from their historical high during the 
period from October 2007 to October 2008. Since November 2008 and until the end of July 2009, the Chinese stock 
markets had been rising dramatically. By calibrating the recent market index price time series to our LPPL model, 
we infer that, in both cases, a bubble had formed in the Chinese stock market and that the market prices were in an 
unsustainable state. We present the an alysis that led us to diagn ose the presence of these two bubbles respectively in 
September 2007 and in July 10, 2009 (iBastiaensen et all 120091) . and to issue an advance notice of the probable time 
of the regime shifts, from a bubble (accelerating "bullish") phase to a ("bearish") regime or a crash. See Fig.Q]for 
an overview of the two bubbles and our predictions. The figure shows the time evolution of two Chinese indexes, the 



'Available at http : //www . csrc . gov . cn/n575458/n4001948/n4002120/4069846 . html accessed on 30 August 2009. 
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dates when we made our predictions and the time intervals of our predicted changes of regime. 

[Figure 1 about here.] 

The organization of the paper is the following. In Sec. 12 we present technical descriptions of all the methods used 
in this paper. Specifically, they are LPPL fitting procedure, Lomb spectral analysis, unit root tests and change-of- 
regime statistic. We present the results of the 2007 and 2009 bubbles in two separate subsections of Sec. [3] In Sec. |4j 
we document and discuss the predictions we made for both bubbles prior to their bursting and, further, describe the 
observations of both markets indicating that these two bubbles actually did burst. Section 5 concludes. 

2. Methods 

Our main method for detecting bubbles and predicting the critical time t c when the bubble will end e ither in a 



crash o r change of reg ime is by fitting observed price time series to a log periodic power law (LPPL) model (I Sornette . 



2003a,b 



Zhou, 



20071) . This is a stochastic fitting procedure that we complement with other techniques, described 
below. This philosophy of using multiple measures aids in filtering predictions, in that a candidate prediction must 
pass all tests to be considered worthy. These techniques form a toolset that has successfully been put to practice over 
the past years by Sornette et al. as described in the introduction. Independently a similar toolset has recently been 
developed within the Research Group of BNP Paribas Fortis (Global Markets) on the same methodology but with a 
slightly different implementation of the fitting procedure and the Lomb analysis. 

2.1. General LPPL fitting technique 

Consider a time series (such as share price) p(t) between starting and ending dates t\ and t%. The LPPL model that 
we use is 

ln[p(0] = A + Bx m + Cx m cos(w In x + <p), ( 1 ) 

where x — t c - 1 measures the time to the critical time t c . For < m < 1 and B < (or m < and B > 0), the power law 
term Bx'" describes the faster-than-exponential acceleration of prices due to positive feedback mechanisms. The term 
proportional to cos(w In x + 4>) expresse s a correction to this super-exponential behavior, which has the symmetry of 



discrete scale invariance ( Sornette, 



19981) . By varying t\ and ?2, we can investigate the stability of the fitting parameters 
with respect to starting and ending points. 

It is worthwhile pointing out that calibrating Eq. (HJ to any given price (or log-price) trajectory will always provide 
some fit parameters. That is, any model can be fit to any data. Hence, it is necessary to establish a constraint — the 
LPPL condition — to filter all of the fitting results. We filter on three parameters: 

t c > h,B < 0,0 < m < 1. (2) 
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This filter selects regimes with faster-than-exponential acceleration of the log-price with a diverging slope at the 
critical future time t c . 

There are four nonlinear parameters (f c , m, to, and cf>) and three linear parameters (A, B, and C) in Eq. (HJ. In order 
to reduce the fitting parameters, the linear parameters are slaved to the nonlinear parameters. By rewriting Eq. (Q]) as 
In p{t) — A + Bf(t) + Cg{t) and using an estimate of the nonlinear parameters, the linear parameters can be solved 
analytically via: 



(3) 
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The i mplementation of the fitting proceeds in two steps. First, we adopt the Taboo search dCviiovic and Klinowski , 



19951) to find 10 candidate solutions from our given search space. Second, each of these solutions is used as an initial 
estimate in a Levenberg-Marquardt nonlinear least squares algorithm. The solution with the minimum sum of squares 
between model and observations is taken as the solution. 



2.2. Stability of fits vs. shrinking and expanding intervals and probabilistic forecasts 

In order to test the sensitivity of variable fitting intervals [/i.fe]. we adopt the strategy of fixing one endpoint and 
varying the other one. For instance, if t% is fixed, the time window shrinks in terms of t\ moving towards ti with a step 
of five days. If t\ is fixed, the time window expands in terms of ti moving away from t\ with a step of five days. For 
each such ?2] interval, the fitting procedure is implemented on the index series three times. Recall that because of 
the rough nonlinear parameter landscape of Eq. ([TJ and the stochastic nature of our initial parameter selection, it is 
expected that each implementation of our fit process will produce a different set of fit parameters. By repeating the 
process multiple times, we investigate an optimal (not necessarily the optimal) region of solution space. 

By sampling many intervals as well as by using bootstrap techniques, we obtain predictions that are inherently 
probabilistic and reflect the intrinsic noisy nature of the underlying generating processes. This allows us to provide 
probabilistic estimations on the time intervals in which a given bubble may end and lead to a new market regime. 
In this respect, we stress that, notwithstanding the common use of the term "crash" to refer to the aftermath of a 
bubble, a real crash does not always occurs. Rather, the end of a bubble may be the most probable time for a crash to 
occur, but the bubble may end without a splash and, instead, transition to a plateau or a slower decay. This point is 
actually crucial in rational expectation models of bubbles in that, even in the presence of investors fully informed of 
the presence of the bubble and with the knowledge of its end da te, it remains rational to stay i nvested in the market to 



garner very large returns since the risk of a crash remains finite (iJohansen et al 



1999 



2000a). 



2.3. Lomb spectral analysis 

Fitting the logarithm of prices to the model Eq. ([]]) gives strong evidence supporting log-periodicity i n that stable 



values of the angular frequency a> are found. We test this feature further by using Lomb spectral analysis ( IPress et al 

5 



199 w for detecting the log-periodic oscillations. The Lomb method is a spectral analysis designed for irregularly sam- 
pled data and gives the same results as the standard Fourier spectral analysis for evenly spaced data. Specifically, given 
a time series, the Lomb analysis returns a series of frequencies a> and associated power at each frequency, Py(fa;)). The 



frequency with the maximum power is taken as the Lomb frequency, WLomb- Following 
spectral Lomb analysis is performed on two types of signals. 

Parametric detrending approach. The first is the series of detrended residuals, calculated as 



Sornette and Zhoul (120021) , the 



r(t) = X - m Qn[p(t)] -A- Bx m ), 



(4) 



where x = t c — t and A, B, t c and m have been found via the method of Section I2TTI ( Joha nsen and SornetteL 11999a; 



Zhou et al 



2007b . As suggested in Eq. 0, the log-periodic oscillations result from the cosine part. The angular 



frequency WLomb is then compared with that found in the LPPL fitting procedure, u^. 

Non-parametric, (H, q) analysis. Th e second is the (H, ap-deriv ative of the logarith mic prices, which has been suc- 



cessfully applied to financial crashes (iZhou and Sornette , 



2003bl) and critical ruptures (IZhou and Sornette t 



the detection of log - period ic components. The (H, q) analysis is a generalization of the g-analysis ( lErzan , 



Erzan and Eckmann 



2002b) for 



1997 



1997b . which is a natural tool for the description of discrete scale invariance. The (H,q)- 
derivative is defined as, 

. f(x) — f(flX) 

(5) 



iqm± m ~ fiqx) 



[(l-q)x] H ■ 

We vary H and q in the ranges [-1, 1] and [0, 1], respectively, and perform the Lomb analysis on the resulting 
series. If H — 1 in Eq. (O, the (H, g)-derivative reduces to the normal ^-derivative, which itself reduces to the 
normal derivative in the limit q — > 1~, Without loss of generality, q is constrained in the open interval (0, 1). The 
advantage of the (H, q) analysis is that there is no need for detrending, as it is automatically accounted for by the finite 
difference and the normalization by the d enominator. This method has been applied for detecting lo g-periodicity 



2002; 



in stock market bubbles and anti-bubbles ( I Sornette a nd Zhou, 
foreign capital inflow bubble end ing in early 2001 (ISornette and Zhou , 



Zhou and Sornette, 



2003b 



2004), in the USA 



bubble (IZhou and Sornette , 



2003a) 



20041) and in the ongoing UK real estate 



2.4. Ornstein- Uhlenbeck and unit root tests 
Recently. 



Lin et al 



(2009) have put forward a self-consistent model for explosive financial bubbles, in which the 
LPPL fitting residuals can be modeled by a mean-reversal Ornstein-Uhlenbeck (O-U) process if the logarithmic price 
in the bubble regime is attributed to a deterministic LPPL component. The test for the O-U property of LPPL fitting 
residuals can be translated into an AR(1) test for the corresponding residuals. Hence, we can verify the O-U property 
of fitting residuals by applying unit-root tests on the residuals. We use the Phillips-Perron and Dickey-Fuller unit-root 



6 



tests. A rejection of null hypothesis Hq suggests that the residuals are stationary and therefore are compatible with the 
O-U process in the residuals. Our tests use the same time windows as the LPPL calibrating procedure. 

2.5. Statistic of change of regime 

We define a simple statistic to demonstrate the change of regime in the "post-mortem" analysis of our prediction. 
We calculate the difference d co (t) between the closing and opening prices on each trading day and track the fraction 
of days with negative d co (t) in a rolling window of T days, with T = 10, 20 and 30 days. 

3. Results 

In the following two subsections, we present our analysis of the separate 2007 and 2009 bubbles using the methods 
described in Section [2] 

3.1. Back test of Chinese bubble from 2005 to 2007 
3.1.1. LPPL fitting with varying window sizes 

As discussed above, we test the stability of fit parameters for the two indexes, SSEC and SZSC, by varying the 
size of the fit intervals. Specifically, the logarithmic index is fit by the LPPL formula, Eq. £[): 

1. in shrinking windows with a fixed end date ?2 = Oct- 10-2007 with the start time t\ increasing from Oct-0 1-2005 
to May-3 1-2007 in steps of five (trading) days and 

2. in expanding windows with a fixed start date t\ = Dec-01-2005 with the end date to increasing from May-01- 
2007 to Oct-01-2007 in steps of five (trading) days. 

In the above two fitting procedures, we fit the indexes 124 times in shrinking windows and 15 times in expanding 
windows. After filtering by the LPPL conditions, we finally observe 72 (78) results in the first step and 11 (15) 
results in the second step for SSEC (respectively SZSC). Figures [2] (a) and (c) illustrate six chosen fitting results of 
the shrinking windows for SSEC and SZSC and Figures [2] (b) and (d) illustrate six fitting results of the expanding 
time intervals for SSEC and SZSC. The dark and light shadow boxes in the figures indicate 20%/80% and 5%/95% 
quantile range of values of the crash dates for the fits that survived filtering. One can observe that the observed market 
peak dates (16 October 2007 for SSEC, 31 October 2007 for SZSC) lie in the quantile ranges of predicted crash dates 
t c using only data from before the market crash (i.e., using t2 < f c .obs)- 

[Figure 2 about here.] 
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3.1.2. Lomb analysis, parametric approach 

Fig.[3]summarizes the results of our Lomb analysis on the detrended residuals r(t). The Lomb periodograms (fV 
with respect to WLomb) are plotted in Fig. [2a) for four typical examples, which are (h,t2) = (Mar- 13-2006, Oct-10- 
2007) and (Dec-12-2005, Sep-07-2007) for SSEC and (Apr-04-2006, Oct-10-2007) and (Dec-01-2005, Sep-09-2007) 
for SZSC. The inset illustrates the corresponding detrended residuals r(t) as a function of ln(f c - t). We select the 
highest peak with its associated WLomb- 

The values of a»Lomb must be consistent with the values of wg t obtained from the fitting. We plot the bivariate 
distribution of pairs («Lomb, P'™ x ) for different LPPL calibrating windows in Fig. [3jb) and find that the minimum 
value of P™ ax is approximately 54 for SSEC and approximately 30 for SZSC. These peaks are linked to a false alarm 
probability, which is defined as the chance that we falsely detect log-periodicity in a signal without true log-periodicity. 

To calculate this false alarm probability, a model of the distribution of the residuals must be used. We 'bracket' 
a range of models, from uncorrec ted white noise to long-range correlated noise. For white noise, we find the false 



alarm probability to be Pr <k 10 (Pr ess et al 



1996). If the residuals have power-law behaviors with exponent in the 
range 2-4 and long-range correlations characterized by a Hurst index H < 0.7, we have Pr < 10 -2 (IZhou and Sornette , 



2002d) . 

The inset of Fig. [3jb) plots Wfi t with respect to d>homb- One can see that most pairs of (wLomb, <^>fu) are overlapping 
on the line y - x, which indicates the consistency between to^ and <x>Lomb- The other pairs are located on the line 
y = 2x. We can interpret these as a fundamental log-periodic component at at and its harmonic compon ent at 2a>. 

periodic components can be expec ted generically in log -periodic signals (ISornette , 



The existence of harmonics of log 



1998 



Gfuzman and Sornette , 



2002 



Zhou and Sornette. 



of financial time series and for other systems (IZhou et al 



2002a, 



2009) and has been documented in earlier studies both 



2007). When the harmonics are well defined with close- 



to-integer ratios to a common fundamental frequency as is the case here, this is in general a diagnostic of a very 
significant log-periodic component. 

[Figure 3 about here.] 

3.1.3. Lomb analysis, non-parametric (H, Q) approach 

In order to non-parametrically check the existence of log-periodicity by means of {H, (^-analysis, we take f(x) = 
\np(t) and x = t c — t with f f = Oct-10-2007 or Oct-25 -2007 (the two observed peak dates of the indexes). For each pair 
of (H, q) values, we calculate the (H, g)-derivative with a given t c , on which we calculate the Lomb analysis. Fig. [3 a) 
illustrates the Lomb periodograms for both indexes with t c = Oct-10-2007, H — 0, q — 0.8 and t c = Oct-25-2007, 
H = 0.5, q = 0.7. The corresponding In p(t), defined by formula (01, is plotted with respect to ln(f e - f) in the 
inset. The highest Lomb peak of the resultant periodogram has height P™ x and abscissa o>Lomb, both P™ ax and WLomb 
being functions of H and q. 

We scan a 21 x 9 rectangular grid in the (H, q) plane, with H e [-1, 1] and q e [0.1,0.9], both in steps of 
0.1. Fig. Htb) illustrates the bivariate distribution of pairs (a>Lomb, ^ > ™ ax )- The inset shows a simple histogram of 



<^Lomb- We observe the three most prominent clusters corresponding to SSEC with t c = Oct-10-2007 (Oct-25-2007) 
as w Lomb = 0.86 + 0.16 (1.06 + 0.41), w Lomb = 4.04 + 0.29 (4.33 + 0.31) and oj 2 Lomh = 10.05 + 0.56 (10.59 ± 0.44) 



For SZSC with t c = Oct-10-2007 (Oct-25-2007), we find the three most prominent clusters to be w° omb = 0.80 + 0.27 
(0.75 + 0.32), Lo l Lomb = 4.21 + 0.38 (4.33 + 0.31), and = 9.29 + 0.29 (9.65 + 0.27). 

The small value of a£ . corresponds to a component with less than one full period within the i nterval of the 



ln(f c — t) variable investigated here. According to extensive tests performed in synthetic time series (iHuang et al 



2000), it should be interpreted as a spurious peak associated with the most probable partial oscillations of a noisy 

l 

Lomb 



signal and/or to a residual global trend in the (H, (^-derivative. Then, col . is identified as the fundamental angular 



log-frequency and <w Lomb ~ 2<w Lomb is interpreted as its second harmonic. 

[Figure 4 about here.] 

3.2. Back test of Chinese bubble from 2008 to 2009 

3.2.1. LPPL fitting with varying window sizes 

The main results of our analysis for the 2008-2009 bubble is illustrated in Fig. [5] The SSEC and SZSC index 
series between Oct-15-2008 and Jul-31-2009 have been calibrated to the LPPL formula given by Eq. ([TJ in shrinking 
and expanding windows. The shrinking windows are obtained by increasing the starting date t\ from Oct-15-2008 to 
Apr-3 1-2009 with a step of five days and keeping the last day ?2 fixed at Jul-31-2009. The expanding windows are 
obtained by fixing the starting day t \ at Nov-01-2008 and moving the ending day to away from t\ from Jun-01-2009 to 
Jul-31-2009 in increments of five days. The results are filtered by the LPPL conditions, resulting in 38 (respectively 
13) fits for SSEC and 28 (respectively 13) fits for SZSC in shrinking windows (respectively increasing windows). 
Figures|5ta) and (c) illustrate six chosen fitting results of the shrinking windows for SSEC and SZSC and Figures[3b) 
and (d) illustrate six fitting results of the increasing time intervals for SSEC and SZSC. The dark and light shadow 
boxes in the figures indicate 20%/80% and 5%/95% quantile range of values of the crash dates for the fits that survived 
filtering. Our calibration confirms the faster- than-exponential growth of the market index over this time interval, a 
clear diagnostic of the presence of a bubble. It also diagnoses that the critical time t c for the end of the bubble and the 
change of market regime lies in the interval August, 1-26, 2009 for SSEC and August, 3-9, 2009 for SZSC (20%/80% 
quantile confidence interval). 

[Figure 5 about here.] 

3.2.2. Lomb analysis, parametric approach 

We use the Lomb spectral analysis technique to further investigate the log-periodic oscillations of Eq. (HJ in both 
indexes from Oct-15-2008 to Jul-31-2009. First, we calculate the detrended residuals r(t) in all the surviving LPPL 
windows and calculate the Lomb periodogram. The highest peak P™ ax and its abscissa «Lomb are extracted from the 
residual Lomb periodograms to plot as points in Fig. HJa). The inset plots Wfi t with respect to WLomb- One can see 
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that most pairs of («Lomb, w fit) are overlapping on the line y — x and the other pairs are located on the line y = 2x, 
confirming the existence of a strong log-periodic component. Fig. [(J a) shows four clusters of values for a>Lomb around 
cj} . = 5 ± 1. oj? . = 10 ± 1, w?~ 4 . = 17.5 ± 2 and , = 27 + 2. The first value a>} . can be interpreted as the 

Lomb ' Lomb ' Lomb Lomb Lomb ' 

main angular log-frequency, while the others are the harmonics of order 2 to 5, with a rather large noise on &»j?~ 4 b anc ' 

^Lomb ' 

3.2.3. Lomb analysis, non-parametric (H, Q) approach 

We next scan a 21 x 9 rectangular grid in the (H, q) plane with H e [-1, 1] and q e [0.1,0.9], both in steps of 
0.1, using t c = Jul-3 1-2009, a date four days before the peak of the SSEC index. We calculate the corresponding 
(H, g)-derivatives In p(t), defined by formula (O, for this set of H and q values. For each obtained In p(t), we 
estimate the Lomb periodogram and plot the highest Lomb peak P™ x as a function of its abscissa WLomb in Fig.|6jb). 
The inset shows the simple histogram of WLomb- 

For the Shanghai index (SSEC), the three most prominent clusters correspond to <^ omb = 1.12 + 0.66, ^ omb = 
4.0 + 1 and cui^ , = 17.6 + 2.7. We interpret the first cluster, <y? as due to the noise decorating the power-law. 
This is because it corresponds to a component with less than one full period within the interval of the l n(t r - 1) variable 



investigated here. According to extensive tests performed in synthetic time series dHuang et al 



peak associated the most probable partial oscillations of a noisy signal. We identify the second cluster, a>} , , as the 



20001) . it is a spurious 

l 

Lomb' 



fundamental angular log-frequency for SSEC. The second cluster around (^ omh is compatible with interpreting it as 
being the third harmonics. It is notable tha t the second ha r monic is not visible in this distribution, a phenomenon 



which has bee n reported for other systems ( Johansen et al 



group analysis (IGluzman and Sornette , 



2002). 



2000b) and can be rationalized from a renormalization 



[Figure 6 about here.] 

3.3. Unit root tests of the 2005-2007 and 2008-2009 Chinese bubbles 

We apply unit root tests to the series of residuals for each [t \ , tj\ interval, where the residuals are calculated by 
subtracting the model from the observations. The goal is to investigate the stationarity of the residuals to determine if 
a mean-reversal Ornstein-Uhlenbeck process is a good model for them. The null hypothesis of the unit-root test is that 
the series being tested is non-stationary. Rejection of the null hypothesis, then, implies that the series (or residuals in 
our case) is stationary. Refer to Section l2~4l for further details. 

We calibrate the LPPL model Eq. (Q} to the SSEC and SZSC indexes in two representative intervals. The first 
interval is from Dec-01-2005 to Oct-10-2007 and the other is from Oct-15-2008 to Jul-31-2009. We scan each interval 
with growing and shrinking windows, as described above, and report the fraction Plppl of these different windows 
that meet the LPPL conditions, Eq. ©. We then calculate the conditional probability that, out of the fraction Plppl of 
windows that satisfy the LPPL condition, the null hypothesis of non-stationarity is rejected for the residuals. Results 
of these tests are shown in Table Q] 
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For the time interval from Dec-01-2005 to Oct-10-2007 for the SSEC index (respectively the SZSC index), the 
fraction Plppl = 56.9% (respectively 67.6%) of the windows satisfy the LPPL conditions. All of the fitting residuals 
of both indexes reject the null hypothesis at significance level 0.01 based on the two tests, implying that the residuals 
are stationary. For the fraction of windows which satisfy the LPPL conditions, the fitting residuals of the SSEC 
(respectively SZSC) index can be regarded as generated by a stationary process at the 99.9% (respectively 99%) 
confidence level. 

For the time interval from Oct- 15-2008 to Jul-3 1-2009, we find that a fraction P LPPL = 94.4% (respectively 74.7%) 
of the windows satisfy the LPPL conditions for SSEC (respectively SZSC). All of the fitting residuals of both indexes 
reject the null hypothesis at significance level 0.001 based on the two tests, implying that the residuals are stationary. 
For the fraction of windows which satisfy the LPPL conditions, the fitting residuals of both indexes can be considered 
as a stationary process at the 99.9% confidence level. 

[Table 1 about here.] 

4. Prior prediction of both crashes 

We make the title of this section pleonastic to emphasize that we predicted both crashes with our techniques before 
the actual dates of the observed peaks in the two indexes. The previous sections have presented more thorough 'post- 
mortem' analyses performed after the observed crashes. This section documents the specific but simpler predictions 
that we announced in advance. 

4.1. 2005-2007 bubble 

Two of us (WXZ and DS) performed a LPPL analysis in early September 2007, which led to (i) a diagnostic of 
an on-going bubble and (ii) the prediction of the end of the bubble in early 2008. One of us (DS) communicated this 
prediction on October 18, 2007 at a prominent hedge-fund conference in Stockholm. The participants, managers of 
top global macro hedge-funds, constitute arguably the best proxy for the academic idealization of "rational investors" 
with access to almost unlimited resources and with the largest existing incentives to motivate themselves to acquire all 
possible relevant information and trade accordingly. These participants responded that the predicted change of regime 
was impossible because, in their opinion, the Chinese government would prevent any turmoil on the Chinese stock 
market until at least the end of the Olympic Games in Beijing (August 2008). After the communication of October 
18, 2007, the Hang Seng China Enterprises Index (HSCEI) reached the historical high 20609.10 on 2 November 
2007. Afterwards, the first valley HSCEI= 15460.72 (-25% from historical high) was reached on 22 Nov 2007 and 
the bottom HSCEI=4792.37 (-77% from historical high) was on 29 Oct 2008. On 19 March 2008, HSCEI=1 1379.91 
was another deep valley. These drops occurred after a six-fold appreciation of the Chinese market from mid-2005 to 
October 2007. 
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4.2. 2008-2009 bubble 

4.2.1. The announcement of the prediction 



2009), in which we 



On 10 July 2009, we submitted our prediction online to the |arXiv . org; dBastiaensen et al. 
gave the 20%/80% (respectively 10%/90%) quantiles of the projected crash dates to be July 17-27, 2009 (respectively 
July 10 - August 10, 2009). This corresponds to a 60% (respectively 80%) probability that the end of the bubble 
occurs and that the change of regime starts in the interval July 17-27, 2009 (respectively July 10 - August 10, 2009). 
Redoing the analysis 5 days later with tj = July 14, 2009, the predictions tightened up with a 80% probability for the 
change of regime to start between July 19 and August 3, 2009 (unpublished). 

[Figure 7 about here.] 

The following paragraph and figure[7]are repr oduced from the online ( and un-refereed) prediction, which is avail- 



able in its initial form at the corresponding URL dBastiaensen et al 



2009) 



The result of the analysis is summarized below in the Figure. We analyzed the Shanghai SSE Com- 
posite Index time series between October 15, 2008 and July 9, 2009. We increased the starting date of 
the LPPL analysis in steps of 15 days while keeping the ending date fixed, resulting in 10 fits. The figure 
shows observations of the SSEC Index as black dots (joined by straight lines) and the LPPL fits as smooth 
lines until the last day of analysis. The y-axis is logarithmically scaled, so that an exponential function 
would appear as a straight line and a power law function with a finite-time singularity would appear with 
a slightly upward curvature. Note that the LPPL fits to the observations exhibit this slightly upward cur- 
vature. The vertical and horizontal dashed lines indicate the date and price of the highest price observed, 
July 6, 2009. Extrapolations of the fits to 100 days beyond July 9, 2009 are shown as lighter dashed 
lines. The darker shaded box with diagonal hatching indicates the 20%/80% quantiles of the projected 
crash dates, July 17-27, 2009. The lighter shaded box with horizontal hatching indicates the range of all 
10 projected crash dates, July 10 - August 10, 2009. These two shaded boxes indicate the most probable 
times (with the associated confidence levels) to expect peak and possible subsequent crash of the Index. 
The parameters of the fit confirm the faster-than- exponential growth of the Shanghai SSE Composite 
Index over this time interval, a clear diagnostic of the presence of a bubble. 

4.2.2. What actually happened 

On July 29, 2009, Chinese stocks suffered their steepest drop since November 2008, with an intraday bottom of 
more than 8% and an open-to-close loss of more than 5%. The market rebounded with a peak on August 4, 2009 
before plummeting the following weeks. The SSEC slumped 22 percent in August, the biggest decline among 89 
benchmark indices tracked world wide by Bloomberg, in stark contrast with being the no. 1 performing index during 
the first half of this year. These striking facts show the detachment of the Chinese equity market from other markets. 
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This bubble was probably nucleated by China's central government's reaction to the global financial crisis. Besides 
announcing the huge stimulus plan on 9 November 2008, a loose monetary policy and regulations caused massive new 
loan issuance as shown in Fig. [8] With overproduction and lower global demands, analysts estimate that up to 50% of 
the increase in credit was used to speculate in equities, property and commodities 0. Rumours of asset bubbles were 
widely heard in the market, but when or if they might crash was unknown as usual. 

Note that the change of regime in the SSEC occurred while the total loan of financial institutions was still growing 
at close to its peak YoY 35% monthly rate. This illustrates that the change of regime has occurred in absence of 
any significant modification of the economic and financial conditions or any visible driving force. This observation, 
which should be surprising to most economists and analysts, is fully expected from the mathematical and statistical 
physics of bifurcations and phase transitions on which our LPPL methodology is based: a possibly vanishingly small 
change of some control parameter may lead to a macroscopic bifurcation or phase transition. Rather than leading to 
an absence of predictability, the accelerating susceptibility of the system associated with the approach towards the 
critical point can be diagnosed, as we have shown. The very clear change of regime documented here provides a 
case-in-point demonstrating this concept of an emergent rupture point characterizing the end of the bubble. 

[Figure 8 about here.] 

Figure |9]presents the evolution with time of the close-open statistic introduced in subsection 12. 5l over the period 
from Jan. 2007 to August 25, 2009. The low (respectively high) values of the index correlate well with the ascending 
(respectively descending) trend of the market. One can also observe the recent remarkably abrupt jump upward of the 
close-open statistic at the time scale T = 10 days, confirming the existence of a sudden change of regime. 

[Figure 9 about here.] 

These events unfolded in a rather bullish atmosphere for the Chinese stock markets. For instance, Bloomberg 
reported on July 30, 2009 that billionaire investor Kenneth Fisher emphasized the great success of China's economy 
compared to the rest of the World and that speculation that the "Chinese government will limit bank loans is un- 
founded." Anecdotal sampling of comments on Chinese online forums suggests a majority of doubters until August 
12, after which a majority endorsed the notion of a change of regime. Several commentators stressed again that pre- 
dictions of Chinese stock markets cannot be correct since China's stock markets are heavily influenced by policies 
(known as a policy market). These comments are similar to the disbelief of hedge-fund managers mentioned in sub- 
section |4T|concerning the prediction of the change of regime at the end of 2007, before the 2008 Beijing Olympics. 

2 see, e.g., BNP Paribas FX Weekly Strategist: China Lending Support (31 July 2009); RBS, Local Markets Asia, Alert China: A savings glut is 
causing problems (3 July 2009) 
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5. Discussion 



We have performed a detailed analysis of two financial bubbles in the Chinese stock markets by calibrating the 
LPPL formula ([TJ to two important Chinese stock indexes, Shanghai (SSEC) and Shenzhen (SZSC) from May 2005 to 
July 2009. Bubbles with the property of faster-than-exponential price increase decorated by logarithmic oscillations 
are observed in two distinct time intervals within the period of investigation for both indexes. The first bubble formed 
in the middle of 2005 and burst in October 2007. The other bubble began in November 2008 and reached a peak in 
early August 2009. 

Our back tests of both bubbles find that the LPPL model describes well the behavior of faster-than-exponential 
increase corrected by logarithmic oscillations in both market indexes. The evidence for the presence of log-periodicity 
is provided by applying Lomb spectral analysis on the detrended residuals and (H, g)-derivative of market indexes. 
Unit-root tests, including the Phillips-Perron test and the Dickey-Fuller test, on the LPPL fitting residuals confirm 
the O-U property and, thus, stationarity in th e residuals, which is in good agreement with the consistent model of 



'explosive' financial bubbles dLin et al 



2009) 



While the present paper presents post-mortem analyses, we 



critical date t c of both bubbles in advance of their demise (ISornette , 



emphasize that we predicted the presence and expected 



2007 



Bastiaensen et al 



2009). These two 



successes prolong the series of favorable outcomes follow ing the prediction of the peak in mid-2004 of the real-estate 



bubbl e in the UK by two of us (|Zhou and SornetteU2003al) . of the peak in mid-2006 of the US housing bubble by two 



of us (IZhou and Sornette . 



2009) 



2006b) and of the peak in July 2008 of the global oil bubble by three of us (Sornette et al 



But not all predictions based on the present methodology have fared so well. In particular, 



Lux 



(2009) and 



Rossei 



(120081) have raised severe objections, following the failure of the we ll-publicized prediction pu blished in 2002 that 



the U.S. stock market would follow a downward log -periodic pattern (ISornette and Zhou . 



2002). How can one make 



sense of these contradictory claims? We summarize the present state of the art as follows. 

1 . Prediction of the end of bubbles should not be confused with predictions based on extrapolation, such 
as those associated with antibubbles. There is a confusion between predicting crashes, on the on e hand, and 



pre dicting the con tinuation of an "antibubble" bearish regime, on the other hand. It seems that both 



and 



Luxld2009l) 



Rosser (2008) amalg amate these two issues, when they focus on the failure of the antibubble prediction in 



Sornette and Zhoul (1 20 02") and conclude that "Sornette and his collaborators failed to forecast future crashes." 



There is indeed a fundamental difference between, on the one hand, (i) the prediction of the end of a bubble 

analyzed here, which is characterized by its critical time t c and, on the other hand, (ii) the extrapolation of an 

"antibubble" pattern. This difference is similar to that between (i) the prediction of the approximate parturition 

time of a foetus on the basis of the recording of key variables obtained during its maturation in the uterus of his 

mother and (ii) the prediction of the death of this individual later in old age from an extrapolation of medical 

variables recorded during his adult life. The former (i) is associated with the maturation phase (the financial 
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bubble versus the uterus-foetus development). It has a rather well-defined critical time which signals the tran- 
sition to a new regime (crash/stagnation/bearish versus birth). In contrast, the latter (ii) may be influenced by 
a variety of factors, particularly exogenous, which may shorten or lengthen the life of the antibubble or of the 
individual. One should thus separate the statistics of successes and failures in the prediction of bubbles on 
the one hand and in the prediction of co ntinuation of LPPL antibu bble patterns on the other hand. This was 



the spirit of the experiment proposed by ISornette and Zhoul (120021) to test the distinct hypothesis that bearish 



regimes following a market peak could be predicted when associated with LPPL patterns. 
2. Intrinsic limits of the prediction of the end of an antibubble. For the reasons just mentioned, it is still an 
open problem to determine when an antibubble ends. Our methods show that one cannot avoid a delay of about 



6 months before identifying the end of an antib ubble (IZho u and Sornette 



for the failure of the 2002 antibubble prediction d Sornette and Zhou , 



2002). 



20051) . This is a partial explanation 



3. Track record of the antibubble method. However, one should not forget that, taken as a distinct class sepa- 
rated from that of diagnosing bubbles and their ends, the predictions based on the antibubb le method can count 



several past successes: (a) on the Nikkei antibubble (IJohansen and Sornette , 



1999b 



nese stock market antibubble ( Zhou and Sornette , 



2000) and (b) on the Chi- 



2004J), in addition to the failure mentioned above. This track 



record is insufficient to conclude. More tests in real time should be performed and rigorous methods developed 
to assess the statistical significance of short cata l ogs of success/failure pre dictions can be applied, base d on 



"roulette" approaches (see Chapter 9 in 



Neyman-Pearson or error diagrams (IMolchan 



Sornette! d2003bl) 



1990 



, Bayes' theorem (IJohansen and Sornette , 



1997). 



2000) and 



4. What we learned from the antibubb le prediction failur e 



that the 2002 antibubble prediction of 



Lux (2009) and 



Rosserl (120081) are right to stress 



Sornette and Zhoul (120021) failed. However, a post-mortem analysis in 



Zhou and Sornettd (120051) has revealed an interesting fact. While the prediction failed w hen the S&P500 is 



value d in U.S. dollars, it becomes quite accurate when expressed in euro or British pounds ( IZhou and Sornette , 



20051) . A plausible interpretation would be that the energetic Fed monetary policy of decreasing its lead rate 
from 6.5% in 2000 to 1 % in 2003 has boosted the stock market in local currency from 2003 on, but has degraded 
the dollar, so that the net effect was that the val ue of the US stock market from an international reference point 



was unfolding as expected from the analysis of 



Sornette and Zhoul (120021) . We do not claim that this changed 



the failure into a success. Instead, it illustrates the effect of monetary feedba cks that have to be include d in 



Zhou and Sornettd (l2006al) . 



improved models incorporating fundamental factors, for instance in the spirit of 
5. Track record for diagnosing bubbles and their ends. Our group has announced advanced prediction (not 
just in retrospect) of bubbles and t heir end ( often a crash). The status of these predictions as of 2002 has been 
discussed in details in Chapter 9 of I Sornettd (|2003b|)'s book. As mentioned a bove, subsequent successes include 



the predictions of the peak in mid-2004 of the re al-estate bubble in the UK b y 



Zhou and Sornette (2003a), of the 



peak in mid-2006 o f the US housing bubb le by IZhou and Sornettd d2006bl) and of the peak in July 2008 of the 



global oil bubble by 



Sornette et al 



(120091) . The present analysis on two bubbles in the Chinese market provide 
15 



additional evidence for the relevance of LPPL patterns in the diagnostic of bubbles. 
6. Recent improvements in methodology. While the core model (or forecasting system) has not changed much 
since the late 1990s, several new developments used here allows us to quantify more accurately both the relia- 
bility and the uncertainties. These improvements include multi-window analysis, probability estimates, and a 
consistent LPPL rational expectation model with mean-reverting residuals. Note that the more recent bubble in 
the Chinese indexes was detected and its end or change of regime was predicted independently by two groups 
(the first four authors from academia on the one hand and the last two authors from industry on the other hand) 
with similar results, showing that the model has been well-documented and can be replicated by industrial 
practitioners. In addition, we stress that the method relies essentially on the competition between the posi- 
tive feedback loop of hig her return anticipations competing with negative feedback spirals of crash expectations 



Ide and Sornette , 



1999) and of 



2002 ), whi c h is at the origin of the acceleration oscillations. In the spirit o flLux and Marchesi 



Gallegati et al 



(2008), this is accounted for in heterogeneous agent mo dels by including nonlin 



ear fundamental investm ent styles competing with non linear momentum trading styles dlde and Sornette , 



2002). 



The initial JLS model of 



Johansen et al 



(11999 



2000al) was based on a conventional neoclassical model assum- 
ing a homogeneous rational agent, but it also enriched thi s set-up by introducing heterogeneous noise traders 
driving a crash hazard rate. More recently,!] 



Lin et al 



(2009) have considered an alternative framework which ex- 



tends the JSL model to account for behavioral herding by using a behavioral stochastic discount factor approach, 
with self-consistent mean-reversal residuals. 

In conclusion, given all the above, we feel this technique is the basis of a prediction platform, which we are 
actively developing, motivated by the conviction that this is the only way to make scientific progress in this deli- 
cate and crucial domain of gr eat societal importance, as illustrated by the 2007-2009 financial and economic crisis 



dSornette and Woodard , 



2009) 
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Figure 1: Evolution of the price trajectories of the SSEC index and the SZSC index over the time interval of this analysis. The solid red lines 
indicate the dates of the respective public announcement of our predictions for the two bubbles (October 18, 2007 and July 10, 2009) while the 
grey zones indicate the 20%/80% confidence intervals for which we forecasted the change of regime. Final closing prices shown in these plots are 
10,614.3 (SZSC) and 2683.72 (SSEC) from September 1, 2009. 
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Figure 2: Daily trajectory of the logarithmic SSEC (a,b) and SZSC (c,d) index from May-01-2005 to Oct-18-2008 (dots) and fits to the LPPL 
formula 0. The dark and light shadow box indicate 20/80% and 5/95% quantile range of values of the crash dates for the fits, respectively. The 
two dashed lines correspond to the minimum date of t\ and the maximum date of t 2 . (a) Examples of fitting to shrinking windows with varied t\ and 
fixed t 2 = Oct-10-2007 for SSEC. The six fitting illustrations are corresponding to h = Sep-30-2005, Dec-05-2005, Feb-13-2006, Apr-24-2006, 
Jan-15-2007, and Mar-12-2007. (b) Examples of fitting to expanding windows with fixed ?i = Dec-01-2005 and varied ti for SSEC. The six fitting 
illustrations are associated with t 2 = Aug-20-2007, Aug-29-2007, Sep-07-2007, Sep-17-2007, Sep-26-2007, Oct-05-2007. (c) Examples of fitting 
to shrinking windows with varied t\ and fixed t 2 = Oct-10-2007 for SZSC. The six fitting illustrations are corresponding to t\ = Sep-30-2005, Dec- 
12-2006, Feb-24-2006, May- 12-2006, Jan-09-2007, and Apr-13-2007. (d) Examples of fitting to expanding windows with fixed t\ = Dec-01-2005 
and varied t 2 for SZSC. The six fitting illustrations are associated with t 2 = Aug-01-2007, Aug-10-2007, Aug-24-2007, Sep-07-2007, Sep-21-2007, 
Oct-08-2007. 
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Figure 3: Lomb tests of the detrending residuals r(t) for SSEC and SZSC. The residuals are obtained from Eq. |4) by substituting different survival 
LPPL calibrating windows with the corresponding fitting results including t c , m, and A. (a) Lomb periodograms for four typical examples, which are 
presented in the legend. The time periods followed the index names represent the LPPL calibrating windows. The inset illustrates the corresponding 
residuals r(t) as a function of ln(? c - f). (b) Bivariate distribution of pairs (wLomb, ^'™ ax ) Ior different LPPL calibrating intervals. Each point in the 
figure stands for the highest peak and its associated angular log-frequency in the Lomb periodogram of a given detrended residual series. The inset 
shows coat as a function of WLomb- 
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Figure 4: Lomb tests of (H, </)-derivative of logarithmic indexes, (a) Lomb periodograms of D^f In p(t) for four typical examples, which are t c = 
Oct-10-2007 with H = and q = 0.8 for SSEC, t c = Oct-25-2007 with H = 0.5 and q = 0.7 for SSEC, t c = Oct-10-2007 with H = and q = 0.8 
for SZSC, and t c = Oct-25-2007 with H = 0.5 and q = 0.7 for SZSC, respectively. The inset shows the corresponding plots of D? In p{t) as a 
function of \n(t c - f). (b) Bivariate distribution of pairs (&*Lomt>, P'™ x ) f° r different pairs of (H, q). Each point corresponds the highest Lomb peak 
and its associated angular log-frequency in the Lomb periodogram of the (H, ^-derivative of logarithmically indexes for a given pair (H, q). The 
inset shows the empirical frequency distribution of WLomb- 
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Figure 5: Daily trajectory of the logarithmic SSEC (a,b) and SZSC (c,d) index from Sep-01-2008 to Jul-31-2009 (dots) and fits to the LPPL formula 
(T}. The dark and light shadow box indicate 20/80% and 5/95% quantile range of values of the crash dates for the fits, respectively. The two dashed 
lines correspond to the minimum date of ti and the fixed date of t 2 ■ (a) Examples of fitting to shrinking windows with varied t \ and fixed f2 = 
Jul-31-2009 for SSEC. The six fitting illustrations are corresponding to f, = Oct-15-2008, Nov-07-2008, Dec-05-2008, Jan-05-2008, Feb-06-2008, 
and Feb-20-2008. (b) Examples of fitting to expanding windows with fixed t\ = Nov-01-2008 and varied t 2 for SSEC. The six fitting illustrations 
are associated with t 2 = Jun-01-2009, Jun-10-2009, Jun-19-2007, Jun-29-2007, Jul-13-2007, Jul-27-2007. (c) Examples of fitting to shrinking 
windows with varied t\ and fixed tj = Jul-31-2009 for SZSC. The six fitting illustrations are corresponding to fi = Oct-15-2008, Nov-03-2008, 
Nov-26-2008, Dec-19-2008, Jan-14-2008, and Jan-23-2008. (d) Examples of fitting to expanding windows with fixed t\ = Dec-01-2005 and varied 
t 2 for SZSC. The six fitting illustrations are associated with t 2 = Jun-01-2009, Jun-10-2009, Jun-19-2007, Jun-29-2007, Jul-13-2007, Jul-27-2007. 
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Figure 6: Detection of log-periodicity in the Chinese bubble from 2008 to 2009. (a) Plots of P™"* with respect to tt>L om b for different LPPL 
calibrating windows. The inset illustrates the dependence of wst on &>Lomb- (b) Bivariate distribution of pairs (<i<Lomb, P™**) for different pairs of 
(H, q) of (H, ^-derivatives In p(t), defined by formula {5}- The inset depicts the empirical frequency distribution of WLomb- 
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Figure 7: Shanghai Composite Index with LPPL result, as presented in the July 10, 2009 arXiv . org submission of Bastiaense n et alj 120091) . 
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Figure 8: (left axis, dots) SSEC compared with Total Loans of Financial Institutions as reported by The People's Bank of China ("Summary of 
Sources & Uses of Funds of Financial Institutions" http://www.pbc.gov.cn/english I (right axis, solid line) YoY % monthly change. This shows 
graphically the widespread belief that the credit growth has fueled the last Chinese equity bubble. 
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Figure 9: Left scale: SSE Composite index from Jan. 2007 to August 25, 2009 (closing price 2980.10). Right scale: fraction of days with negative 
(close-open) in moving windows of length T = 10 days (continuous blue line), T = 20 days (dashed green line) and T = 30 days (dotted-dashed 
red line). 
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Table 1: Unit-root tests on the LPPL fitting residuals for SSEC and SZSC index in our two calibrating ranges. Plppl denotes the fraction of 
windows that satisfy the LPPL condition. P stationaryResi.|LPPL denotes the conditional probability that, out of the fraction Plppl of windows that 
satisfy the LPPL condition, the null unit test for non-stationarity is rejected for the residuals. 



index 


calibrating 


number of 


^LPPL 


signif. 


percentage 


of rejecting 


^StationaryResi. |LPPL 


range 


windows 


level 


Phillips-Perron 


Dickery-Fuller 


SSEC 


2005/12/01 


146 


56.9% 


a = 0.01 


100.0% 


100.0% 


100.0% 


2007/10/10 


a = 0.001 


95.2% 


95.2% 


100.0% 


SZSC 


2005/12/01 


139 


67.6% 


a = 0.01 


100.0% 


100.0% 


100.0% 


2007/10/10 


a = 0.001 


81.3% 


81.3% 


72.3% 


SSEC 


2008/10/15 


54 


94.4% 


a = 0.01 


100.0% 


100.0% 


100.0% 


2009/07/31 


a = 0.001 


100.0% 


100.0% 


100.0% 


SZSC 


2008/10/15 


54 


74.7% 


a = 0.01 


100.0% 


100.0% 


100.0% 


2009/07/31 


a = 0.001 


100.0% 


100.0% 


100.0% 
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